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Abstract 

Mexico is one of the most volcanically active regions 
in North America. Volcanic activity in central Mex- 
ico is associated with the subduction of the Cocos 
and Rivera plates beneath the North American plate. 
Periods of enhanced microseismic activity, associated 
with the volcanic activity of the Popocatepetl vol- 
cano is compared with periods, during which the mi- 
croseismic activity was low. We detected systemati- 
cal changes in the number of lineaments, associated 
with the microseismic activity due to lineament anal- 
ysis of a temporal sequence of high resolution satellite 
images of the Popocatepetl volcano, provided by the 
ASTER/VNIR instrument. The Lineament Extrac- 
tion and Stripes Statistic Analysis (LESSA) software 
package was used for the lineament extraction. In 
the future it would allow develop a methodology for 
detection of possible elevation of pressure in volcano 
edifice. 



1 Introduction 

Popocatepetl is the most famous volcano in Mexico 
located about 70 km to the southeast of Mexico City 
(Long. 261.37, Lat. 19.07, elevation 5465 m). It is 
one of several active volcanoes that form the Trans- 



Volcanic Belt of Mexico (also know as Neo- Volcanic 
Axes) and its existence is related to the geodynamics 
of the North American and Coco plates. It has an 
almost symmetric cone that rises 4,200 m above the 
surrounding plain. The volcano is composed of two 
distinct structures: the old part, called Nexpayantla, 
and the newer structure, which stands on the sout- 
hearstern part of the old Popocatepetl. The crater 
is 400 by 600 m wide with vertical walls. Preva- 
lent volcanic activities, observed in the volcano, are 
Plinian, Vulcanian, pheatomagmatic, lava flows and 
fumaroles. The most recent eruption began in De- 
cember 1994. 

Similarly to other volcanoes, the Popocatepetl vol- 
cano seismicity is represented by of micro- or small 
earthquakes in form of swarms, the magnitude of 
which is generally below 4.0. Presence of earthquake 
swarms is normally related to volcanic activity even if 
the surface eruptions are not observed. Moreover, the 
migration of earthquakes usually suggests the poten- 
tial movement of magma toward the surface. How- 
ever, only a few studies have been made to establish 
a relationship between the appearance and migration 
of hypocenters and the volcanic processes of intru- 
sions. In particular Gambino et al., (2004) analyzed 
the location of microseismicity preceding the 2002- 
2003 Mt. Etna eruption. They found that seismic- 
ity and swarms occurring since February 2002 near 



the central crater can be interpreted as a response of 
the volcano to magma movements along the NNW- 
SSE volcano-genetic trend. In particular, the intense 
micro-seismic activity observed between August 31 
and September 3, 2002 beneath the SSE Rift and 
eastward from it may be related to rising of magma. 

Earthquake swarms and ground deformation arc 
amongst the most frequently cited precursors of vol- 
canic eruptions (Scarpa and Tilling, 1996). In par- 
ticular, Voight (1988, 1989) suggested a forecast- 
ing model based on dynamics of material failure 
to describe accelerating earthquake event rate and 
groimd deformation. It is based on the assumption 
that the growth of magmatic pathways is driven by 
rock failure in the volcanic edifice under sustained, 
near-constant fluid pressure in the magma chamber. 
Therefore the rate of magma ascent is limited by 
the rate of fracture growth, recorded as seismicity. 
Recently Chastin and Main (2003) have tested the 
Voigts hypothesis using the seismic data from vol- 
cano Kilauea, Hawaii. They have found that ac- 
celeration of earthquake rates is not observed sys- 
tematically before all individual eruptions. It can 
be related to the fact that there are a number of 
processes leading to the accumulation of strain and 
therefore to the ground deformation. When the ini- 
tial period of volcano-wide deformation is generally 
attributed to the inflation of the volcanic ediflce by 
the rising magma, the thermal expansion and fault- 
ing/fracturing can also contribute to the strain accu- 
mulation (expansion) o relaxation (fracturing) (Jack- 
son et al., 1998). 

During the last decade, the application of Syn- 
thetic Aperture Radar (SAR) interfcrometry data for 
ground deformation studies on volcanoes allowed to 
generate deformation time series which significantly 
contribute to better understanding of its sources 
(see, for example, Schmidt and Burgmann, (2003), 
Pritchard and Simons, (2002)). Remy et al (2007) 
studied the inflation of the Aira Caldera (Japan) us- 
ing the SAR interfcrometry ERS data during a three 
year period. They estimate a roughly constant an- 
nual rate of inflation of about m3/year that corre- 
sponds to the ground deformation of about 2 cm, 
which agrees with JERS and GPS measurements. 

SAR interfcrometry technique has also been suc- 



cessfully applied to study of ground deformations as- 
sociated with earthquakes related to the movement 
of tectonic plates (see for example (Satybala, 2006; 
Schmidt and Brgmann, 2006, Lasserrc et al., 2005, 
Funning et al., 2005)). Observed deformations were 
of the same order of magnitude as in case of vol- 
cano deformation. Recently Arellano ct al. (2004, 
2006, 2007) have found that presence of earthquake 
can strongly affect the system of lineaments extracted 
from the high- resolution ASTER satellite images. It 
was found that a significant number of lineaments ap- 
pears approximately one month before strong earth- 
quake. One month after the earthquake the linea- 
ment configuration returns to its initial state. These 
features were detected during the analysis of 8 strong 
earthquakes in South and North America and China. 
These features were not observed in tests areas, situ- 
ated far from any earthquake. Singh V.P and Singh 
R.P. (2005) used the lineament analysis to study 
changes in stress pattern around the epicenter of the 
Mw=7.6 Bhuj earthquake occurred 26 January 2001 
in India. The results obtained also confirm that the 
lineaments retrieved from the images 22 days before 
the earthquake differ from the lineaments obtained 3 
days after the earthquake. It was assumed that they 
are related to fractures and faults and that their ori- 
entation and density give an idea about the fracture 
pattern of rocks. The results also show the high level 
of correlation between the continued horizontal max- 
imum compressive stress deduced from the lineament 
and the earthquake focal mechanism. 

The main question is how the lineaments extracted 
from 15-30 m resolution images are able to reflect the 
accumulation of stress deep in the earth crust when 
the ground deformation, associated with earthquakes 
is about a few centimeters? The term lineaments, 
originally proposed by Hobbs (1904), defined linea- 
ments as significant lines of landscape that reveal the 
hidden architecture of the rock basement. However 
O'Leary et. al. (1976) defined lineament as a sim- 
ple or composite linear feature of a surface whose 
parts are aligned in a rectilinear or slightly curvilin- 
ear relationship and which differ from the pattern of 
adjacent features and presumably refiects some sub- 
surface phenomenon. It supposes that lineaments are 
able to detect, at least partially, the presence of rup- 
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tures deep in the Earths crust. This means that cor- 
responding image processing of extraction of system 
of lineaments using different algorithms (see, for ex- 
ample, Karnieli et al., 1996, Fitton and Fox, 1998, 
Wang et al., 1990, Zlatopolsky, 1992, 1997) makes 
possible to retrieve (at least partially) the informa- 
tion about the dynamics of the Earths crust kilome- 
ters and even tens of kilometers underground that is 
disseminated over a surface. Deeper events are dis- 
seminated over bigger areas and require lower reso- 
lution images for their detection. That is why 10-30 
m resolution images are useful for earthquake stud- 
ies. Although they are not able to see some individ- 
ual crack, they are able to integrate (at least) the 
information about the presence of faults tens of kilo- 
meters underground, and track the changes related 
to accumulation or relaxation of strength due to the 
movement of tectonic plates. 

In this paper we will use the same methodology 
to detect possible changes in the lineament structure 
related to the microseismicity events related to the 
volcanic activity of the Popocatepetl volcano. 



2 Instrumentation and Data 
Analysis 

For this study we used the images from the Ad- 
vanced Spaceborne Thermal Emission and Reflection 
Radiometer (ASTER) on board of the TERRA satel- 
lite. The satellite was launched to a circular solar- 
synchronous orbit with altitude of 705 km. The ra- 
diometer is composed by three instruments: Visible 
and Near Infrared Radiometer (VNIR) with 15 m res- 
olution (bands 1-3), Short Wave Infrared Radiometer 
(SWIR) with 30 m resolution (bands 4-9) and Ther- 
mal Infrared Radiometer TIR with 90 m resolution 
(bands 11-14) which measure the reflected and emit- 
ted radiation of the Earths surface covering the range 
0.56 to 11.3 ^nJl (Abrams, 2000). 

The images were processed using the Lineament 
Extraction and Stripes Statistic Analysis (LESSA) 
software package (Zlatopolsky, 1992, 1997), which 
provides a statistical description of the position and 
orientation of short linear structures through detec- 



tion of small linear features (stripes) and calculation 
of descriptors that characterize the spatial distribu- 
tion of stripes. The program also makes it possible 
to extract the lineaments straight lines crossing a sig- 
nificant part of the image. To make this extraction, 
a set of very long and very narrow (a few pixels) 
windows (bands), crossing the entire image in differ- 
ent directions, was used. In each band the density 
of stripes, the direction of which is coincident with 
the direction of the band, is calculated. When the 
density of stripes overcomes a pre-established thresh- 
old, the chain of stripes along the band is considered 
as a lineament. The value of threshold depends on 
the brightness of the image, relief, etc. and is es- 
tablished empirically. Previous studies showed that 
lineaments, extracted from the image by applying the 
LESSA program, are strongly related to the main lin- 
eaments, obtained from the geomorphological studies 
(Zlatopolsky, 1992, 1997). The details about the ap- 
plication of LESSA package for earthquake studies 
are given in (Arellano et al., 2006). 

We studied a time sequence of images covering the 
Popocatepetl volcano area during 2000-2005. Unfor- 
tunately the level of cloudiness was high, and from 10 
images analysed we selected only four of them cover- 
ing exactly the same area with percents of clouds 
(see Figure 1 first line). The area extracted from 
the VNIR band ILevel 2 images correspond to the 
north-east part of the volcano and have the coordi- 
nates, upper left corner: 19°07'56"N, 98°37'20"W; 
upper right corner: 19°07'56"N, 983007W; lower 
right corner: 19°00'55"N, 98°30'07"W; lower left 
corner:19°00'55"N, 98°37'20"W. The corresponding 
systems of lineaments were extracted by applying the 
LESSA program, using the same parameters (Figure 
1, second line). As it can be seen, there are sim- 
ilarities in the structure of lineaments, nevertheless 
their density changes significantly from image to im- 
age, such that we ordered the image sequence from 
high to low lineament density. Subsequent lines con- 
tain the information about the depth, longitude, and 
latitude of microseismic events 30 days before (grey) 
and 90 days after (black) of the time of each image. 
The last line represents the three-dimensional view 
of the same events. These data were provided by 



CENAPRED (http://www.cenapred.unam.mx/es/ 1 
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Figure 1: From top to bottom: fragments of ASTER/VNIR image of the Popocatepetl volcano (band 1), 
systems of lineaments, extracted from the images above, depth, latitude, and longitude of micro earthquakes 
30 davs before fgrav) and 90 davs after (black) to the time of each image, three-dimensional distribution of 
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Figure 2: Relationship between the density of hn- 
eaments and the number of earthquakes during next 
10, 20, 30, 40 and 50 days. 

As it can be seen, there is a relationship between 
the lineament density and the number of earthquakes. 
In order to quantify this relationship at least in very 
elemental way, we plot a relationship between the 
percentage of pixels corresponding to the lineaments 
and the number of earthquakes during the next 10, 
20, 30, 40 and 50 days (Figure 2). To our surprise, a 
very interesting inverse power low dependence with 
an exponent -2.030.07 appeared between the percent- 
age of pixels and the number of earthquakes ten days 
ahead. After that this relationship starts to saturate 
gradually. This is rather surprising and has to the 
best of our never been described. We do not have 
at this point any explanation to offer leaving these 
results as a prospective goal to be investigated. 

3 Discussion and conclusions 

Analysis of a relationship between the structure 
of lineaments and microseismicity observed in the 
Popocatepetl volcano, showed a power low relation- 
ship between the percentage of pixels corresponding 
of lineaments extracted from the ASTER satellite im- 
age and the number of earthquakes ten days ahead. 



This result is opposite to that obtained previously 
for earthquakes not related to the volcano activity. 
In that case the number of lineaments increased sig- 
nificantly a few weeks before an earthquake (Arellano 
et al. 2004, 2006, 2007). This discrepancy can be ex- 
plained assuming that in the last case the main rea- 
son of earthquakes is compression and accumulation 
of strength in the Earths crust due to subduction of 
tectonic plates, whereas in the first case we deal with 
the infiation of a volcano edifice due to elevation of 
pressure and magma intrusion. So it is natural to 
expect that the lineament behavior would be differ- 
ent. This study represents a first step in understand- 
ing of how the processes of volcanic seismic activity 
are reflected in the lineament structure. It is nec- 
essary to increase significantly the number of events 
analysed and develop more sophisticated methodol- 
ogy including both more precise quantification of lin- 
eament behavior, and more precise analysis of seis- 
mic events, including their location, energy release 
and waveform. Nevertheless, the results obtained are 
relevant because open the possibility of volcano moni- 
toring using satellite images. It is important because 
to date only a very low percentage of volcanoes is 
constantly monitored using very efficient but expen- 
sive ground chains. The use of new techniques would 
allow in the future to detect possible elevation of pres- 
sure in the volcano edifice that makes it possible to 
concentrate on a special attention to this specific vol- 
cano and to initiate other studies using more precise 
but unfortunately much more expensive satellite and 
ground based techniques. 
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